{smcl}
{com}{sf}{ul off}{txt}{.-}
      name:  {res}<unnamed>
       {txt}log:  {res}N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\log_Table5.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res} 6 Jul 2024, 21:00:23
{txt}
{com}. 
. //Running the Negative Binomial Regresions (Table 5 in the paper)
. use "N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\complete_dataset.dta", clear
{txt}
{com}. 
. nbreg catch_sf ssb_sf c.bt_tmp##c.bt_tmp  if (draw==100 & state=="NJ"), vce(robust)

{txt}Fitting Poisson model:
{res}
{txt}Iteration 0:{space 2}Log pseudolikelihood = {res:  -4975285}  
Iteration 1:{space 2}Log pseudolikelihood = {res:  -4975285}  
{res}
{txt}Fitting constant-only model:
{res}
{txt}Iteration 0:{space 2}Log pseudolikelihood = {res:-2970710.7}  
Iteration 1:{space 2}Log pseudolikelihood = {res:-2945628.1}  
Iteration 2:{space 2}Log pseudolikelihood = {res:-2945628.1}  
{res}
{txt}Fitting full model:
{res}
{txt}Iteration 0:{space 2}Log pseudolikelihood = {res:-2942728.2}  
Iteration 1:{space 2}Log pseudolikelihood = {res:-2942706.1}  
Iteration 2:{space 2}Log pseudolikelihood = {res:-2942706.1}  
{res}
{txt}{col 1}Negative binomial regression{col 54}{lalign 13:Number of obs}{col 67} = {res}{ralign 9:1,146,957}
{txt}{col 54}{lalign 13:Wald chi2({res:3})}{col 67} = {res}{ralign 9:5218.37}
{txt}{col 1}Dispersion: {res:mean}{col 54}{lalign 13:Prob > chi2}{col 67} = {res}{ralign 9:0.0000}
{txt}{col 1}{lalign 20:Log pseudolikelihood}{col 21} = {res}{ralign 10:-2942706.1}{txt}{col 54}{lalign 13:Pseudo R2}{col 67} = {res}{ralign 9:0.0010}

{txt}{hline 18}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 19}{c |}{col 31}    Robust
{col 1}         catch_sf{col 19}{c |} Coefficient{col 31}  std. err.{col 43}      z{col 51}   P>|z|{col 59}     [95% con{col 72}f. interval]
{hline 18}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 11}ssb_sf {c |}{col 19}{res}{space 2} .1521591{col 31}{space 2} .0032834{col 42}{space 1}   46.34{col 51}{space 3}0.000{col 59}{space 4} .1457237{col 72}{space 3} .1585945
{txt}{space 11}bt_tmp {c |}{col 19}{res}{space 2} .2144025{col 31}{space 2} .0090221{col 42}{space 1}   23.76{col 51}{space 3}0.000{col 59}{space 4} .1967196{col 72}{space 3} .2320855
{txt}{space 17} {c |}
c.bt_tmp#c.bt_tmp {c |}{col 19}{res}{space 2}-.0140455{col 31}{space 2} .0004895{col 42}{space 1}  -28.70{col 51}{space 3}0.000{col 59}{space 4}-.0150048{col 72}{space 3}-.0130862
{txt}{space 17} {c |}
{space 12}_cons {c |}{col 19}{res}{space 2} .5616475{col 31}{space 2} .0405132{col 42}{space 1}   13.86{col 51}{space 3}0.000{col 59}{space 4} .4822432{col 72}{space 3} .6410519
{txt}{hline 18}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}/lnalpha {c |}{col 19}{res}{space 2} .3580541{col 31}{space 2} .0017023{col 59}{space 4} .3547176{col 72}{space 3} .3613905
{txt}{hline 18}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
            alpha {c |}{col 19}{res}{space 2} 1.430543{col 31}{space 2} .0024352{col 59}{space 4} 1.425778{col 72}{space 3} 1.435324
{txt}{hline 18}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. etable,  title("Negative Binomial Reg. Model (NJ)")  note("Robust Standard Errors") notestyles(font(Calibri, size(6))) mstat(N) mstat(r2_p, nformat(%5.4f) label("Pseudo R2")) export("N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\Table5_NJ.txt", replace)
{res}
{smcl}
{reset}{...}
{p}Negative Binomial Reg. Model (NJ){p_end}
{hline 22}{c -}{hline 8}
{space 22} catch_sf
{hline 22}{c -}{hline 8}
ssb_sf{space 16} {space 3}{result:0.152}
{space 22} {space 1}{result:(0.003)}
bt_tmp{space 16} {space 3}{result:0.214}
{space 22} {space 1}{result:(0.009)}
bt_tmp # bt_tmp{space 7} {space 2}{result:-0.014}
{space 22} {space 1}{result:(0.000)}
Intercept{space 13} {space 3}{result:0.562}
{space 22} {space 1}{result:(0.041)}
lnalpha{space 15} {space 3}{result:0.358}
{space 22} {space 1}{result:(0.002)}
alpha{space 17} {space 3}{result:1.431}
{space 22} {space 1}{result:(0.002)}
Number of observations {space 1}{result:1146957}
Pseudo R2{space 13} {space 2}{result:0.0010}
{hline 22}{c -}{hline 8}
{p}Robust Standard Errors{p_end}
{res}{txt}{p 0 1 2}
(collection {res:ETable} exported to file {browse "N:/Faculty/JHolzer/ado/McConnell_Holzer_RDUT/Negative_binomial_regression/Table5_NJ.txt":N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\Table5_NJ.txt})
{p_end}

{com}. putexcel set "N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\nbin_param_nj.xlsx", replace
{res}{p}{txt}note: file will be replaced when the first {cmd:putexcel} command is issued.{p_end}

{com}.     putexcel A1 = matrix(e(b))
{res}{txt}file {bf:N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\nbin_param_nj.xlsx} saved

{com}.         
. nbreg catch_sf ssb_sf c.bt_tmp##c.bt_tmp  if (draw==100 & state=="NY"), vce(robust)

{txt}Fitting Poisson model:
{res}
{txt}Iteration 0:{space 2}Log pseudolikelihood = {res:-2836231.3}  
Iteration 1:{space 2}Log pseudolikelihood = {res:-2836231.3}  
{res}
{txt}Fitting constant-only model:
{res}
{txt}Iteration 0:{space 2}Log pseudolikelihood = {res:-1852406.3}  
Iteration 1:{space 2}Log pseudolikelihood = {res:-1843053.4}  
Iteration 2:{space 2}Log pseudolikelihood = {res:-1843051.4}  
Iteration 3:{space 2}Log pseudolikelihood = {res:-1843051.4}  
{res}
{txt}Fitting full model:
{res}
{txt}Iteration 0:{space 2}Log pseudolikelihood = {res:-1840097.3}  
Iteration 1:{space 2}Log pseudolikelihood = {res:-1840070.7}  
Iteration 2:{space 2}Log pseudolikelihood = {res:-1840070.7}  
{res}
{txt}{col 1}Negative binomial regression{col 56}{lalign 13:Number of obs}{col 69} = {res}{ralign 7:760,485}
{txt}{col 56}{lalign 13:Wald chi2({res:3})}{col 69} = {res}{ralign 7:5177.72}
{txt}{col 1}Dispersion: {res:mean}{col 56}{lalign 13:Prob > chi2}{col 69} = {res}{ralign 7:0.0000}
{txt}{col 1}{lalign 20:Log pseudolikelihood}{col 21} = {res}{ralign 10:-1840070.7}{txt}{col 56}{lalign 13:Pseudo R2}{col 69} = {res}{ralign 7:0.0016}

{txt}{hline 18}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 19}{c |}{col 31}    Robust
{col 1}         catch_sf{col 19}{c |} Coefficient{col 31}  std. err.{col 43}      z{col 51}   P>|z|{col 59}     [95% con{col 72}f. interval]
{hline 18}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 11}ssb_sf {c |}{col 19}{res}{space 2} .2557659{col 31}{space 2} .0043872{col 42}{space 1}   58.30{col 51}{space 3}0.000{col 59}{space 4} .2471671{col 72}{space 3} .2643646
{txt}{space 11}bt_tmp {c |}{col 19}{res}{space 2} .1569635{col 31}{space 2} .0181043{col 42}{space 1}    8.67{col 51}{space 3}0.000{col 59}{space 4} .1214799{col 72}{space 3} .1924472
{txt}{space 17} {c |}
c.bt_tmp#c.bt_tmp {c |}{col 19}{res}{space 2}-.0115206{col 31}{space 2} .0009796{col 42}{space 1}  -11.76{col 51}{space 3}0.000{col 59}{space 4}-.0134406{col 72}{space 3}-.0096006
{txt}{space 17} {c |}
{space 12}_cons {c |}{col 19}{res}{space 2} .5777214{col 31}{space 2} .0826598{col 42}{space 1}    6.99{col 51}{space 3}0.000{col 59}{space 4} .4157112{col 72}{space 3} .7397317
{txt}{hline 18}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}/lnalpha {c |}{col 19}{res}{space 2} .2796915{col 31}{space 2}  .002387{col 59}{space 4} .2750131{col 72}{space 3} .2843699
{txt}{hline 18}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
            alpha {c |}{col 19}{res}{space 2} 1.322722{col 31}{space 2} .0031573{col 59}{space 4} 1.316548{col 72}{space 3} 1.328924
{txt}{hline 18}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. etable,  title("Negative Binomial Reg. Model (NY)")  note("Robust Standard Errors") notestyles(font(Calibri, size(6))) mstat(N) mstat(r2_p, nformat(%5.4f) label("Pseudo R2")) export("N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\Table5_NY.txt", replace)
{res}
{smcl}
{reset}{...}
{p}Negative Binomial Reg. Model (NY){p_end}
{hline 22}{c -}{hline 8}
{space 22} catch_sf
{hline 22}{c -}{hline 8}
ssb_sf{space 16} {space 3}{result:0.256}
{space 22} {space 1}{result:(0.004)}
bt_tmp{space 16} {space 3}{result:0.157}
{space 22} {space 1}{result:(0.018)}
bt_tmp # bt_tmp{space 7} {space 2}{result:-0.012}
{space 22} {space 1}{result:(0.001)}
Intercept{space 13} {space 3}{result:0.578}
{space 22} {space 1}{result:(0.083)}
lnalpha{space 15} {space 3}{result:0.280}
{space 22} {space 1}{result:(0.002)}
alpha{space 17} {space 3}{result:1.323}
{space 22} {space 1}{result:(0.003)}
Number of observations {space 2}{result:760485}
Pseudo R2{space 13} {space 2}{result:0.0016}
{hline 22}{c -}{hline 8}
{p}Robust Standard Errors{p_end}
{res}{txt}{p 0 1 2}
(collection {res:ETable} exported to file {browse "N:/Faculty/JHolzer/ado/McConnell_Holzer_RDUT/Negative_binomial_regression/Table5_NY.txt":N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\Table5_NY.txt})
{p_end}

{com}. putexcel set "N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\nbin_param_ny.xlsx", replace
{res}{p}{txt}note: file will be replaced when the first {cmd:putexcel} command is issued.{p_end}

{com}.     putexcel A1 = matrix(e(b))
{res}{txt}file {bf:N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\nbin_param_ny.xlsx} saved

{com}.         
. nbreg catch_sf ssb_sf c.bt_tmp##c.bt_tmp  if (draw==100 & state=="CT"), vce(robust)

{txt}Fitting Poisson model:
{res}
{txt}Iteration 0:{space 2}Log pseudolikelihood = {res:-377943.05}  
Iteration 1:{space 2}Log pseudolikelihood = {res:-377943.04}  
{res}
{txt}Fitting constant-only model:
{res}
{txt}Iteration 0:{space 2}Log pseudolikelihood = {res:-256648.19}  
Iteration 1:{space 2}Log pseudolikelihood = {res:-254569.56}  
Iteration 2:{space 2}Log pseudolikelihood = {res: -254565.8}  
Iteration 3:{space 2}Log pseudolikelihood = {res: -254565.8}  
{res}
{txt}Fitting full model:
{res}
{txt}Iteration 0:{space 2}Log pseudolikelihood = {res:-254067.15}  
Iteration 1:{space 2}Log pseudolikelihood = {res:-254062.47}  
Iteration 2:{space 2}Log pseudolikelihood = {res:-254062.47}  
{res}
{txt}{col 1}Negative binomial regression{col 56}{lalign 13:Number of obs}{col 69} = {res}{ralign 7:117,865}
{txt}{col 56}{lalign 13:Wald chi2({res:3})}{col 69} = {res}{ralign 7:965.61}
{txt}{col 1}Dispersion: {res:mean}{col 56}{lalign 13:Prob > chi2}{col 69} = {res}{ralign 7:0.0000}
{txt}{col 1}{lalign 20:Log pseudolikelihood}{col 21} = {res}{ralign 10:-254062.47}{txt}{col 56}{lalign 13:Pseudo R2}{col 69} = {res}{ralign 7:0.0020}

{txt}{hline 18}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 19}{c |}{col 31}    Robust
{col 1}         catch_sf{col 19}{c |} Coefficient{col 31}  std. err.{col 43}      z{col 51}   P>|z|{col 59}     [95% con{col 72}f. interval]
{hline 18}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 11}ssb_sf {c |}{col 19}{res}{space 2} .1775846{col 31}{space 2} .0137416{col 42}{space 1}   12.92{col 51}{space 3}0.000{col 59}{space 4} .1506515{col 72}{space 3} .2045177
{txt}{space 11}bt_tmp {c |}{col 19}{res}{space 2} .6221473{col 31}{space 2} .0304483{col 42}{space 1}   20.43{col 51}{space 3}0.000{col 59}{space 4} .5624696{col 72}{space 3} .6818249
{txt}{space 17} {c |}
c.bt_tmp#c.bt_tmp {c |}{col 19}{res}{space 2}-.0327126{col 31}{space 2} .0015141{col 42}{space 1}  -21.60{col 51}{space 3}0.000{col 59}{space 4}-.0356803{col 72}{space 3}-.0297449
{txt}{space 17} {c |}
{space 12}_cons {c |}{col 19}{res}{space 2}-2.061887{col 31}{space 2} .1486521{col 42}{space 1}  -13.87{col 51}{space 3}0.000{col 59}{space 4} -2.35324{col 72}{space 3}-1.770535
{txt}{hline 18}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}/lnalpha {c |}{col 19}{res}{space 2} .3568311{col 31}{space 2} .0064893{col 59}{space 4} .3441123{col 72}{space 3} .3695498
{txt}{hline 18}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
            alpha {c |}{col 19}{res}{space 2} 1.428794{col 31}{space 2} .0092718{col 59}{space 4} 1.410737{col 72}{space 3} 1.447083
{txt}{hline 18}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. etable,  title("Negative Binomial Reg. Model (CT)")  note("Robust Standard Errors") notestyles(font(Calibri, size(6))) mstat(N) mstat(r2_p, nformat(%5.4f) label("Pseudo R2")) export("N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\Table5_CT.txt", replace)
{res}
{smcl}
{reset}{...}
{p}Negative Binomial Reg. Model (CT){p_end}
{hline 22}{c -}{hline 8}
{space 22} catch_sf
{hline 22}{c -}{hline 8}
ssb_sf{space 16} {space 3}{result:0.178}
{space 22} {space 1}{result:(0.014)}
bt_tmp{space 16} {space 3}{result:0.622}
{space 22} {space 1}{result:(0.030)}
bt_tmp # bt_tmp{space 7} {space 2}{result:-0.033}
{space 22} {space 1}{result:(0.002)}
Intercept{space 13} {space 2}{result:-2.062}
{space 22} {space 1}{result:(0.149)}
lnalpha{space 15} {space 3}{result:0.357}
{space 22} {space 1}{result:(0.006)}
alpha{space 17} {space 3}{result:1.429}
{space 22} {space 1}{result:(0.009)}
Number of observations {space 2}{result:117865}
Pseudo R2{space 13} {space 2}{result:0.0020}
{hline 22}{c -}{hline 8}
{p}Robust Standard Errors{p_end}
{res}{txt}{p 0 1 2}
(collection {res:ETable} exported to file {browse "N:/Faculty/JHolzer/ado/McConnell_Holzer_RDUT/Negative_binomial_regression/Table5_CT.txt":N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\Table5_CT.txt})
{p_end}

{com}. 
. putexcel set "N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\nbin_param_ct.xlsx", replace
{res}{p}{txt}note: file will be replaced when the first {cmd:putexcel} command is issued.{p_end}

{com}.     putexcel A1 = matrix(e(b))
{res}{txt}file {bf:N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\nbin_param_ct.xlsx} saved

{com}. 
. log close
      {txt}name:  {res}<unnamed>
       {txt}log:  {res}N:\Faculty\JHolzer\ado\McConnell_Holzer_RDUT\Negative_binomial_regression\log_Table5.smcl
  {txt}log type:  {res}smcl
 {txt}closed on:  {res} 6 Jul 2024, 21:36:46
{txt}{.-}
{smcl}
{txt}{sf}{ul off}